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Using the recently found by G. Horowitz and M. Roberts (arXiv:0908.3677) numerical model of 
the ground state of holographic superconductors (at zero temperature), we calculate the conductiv- 
ity for such models. The universal relation connecting conductivity with the reflection coefficient 
was used for finding the conductivity by the WKB approach. The dependence of the conductivity 
on the frequency and charge density is discussed. Numerical calculations confirm the general argu- 
ments of (arXiv:0908.3677 ) in favor of non-zero conductivity even at zero temperature. In addition 
to the Horowitz-Roberts solution we have found (probably infinite) set of extra solutions which are 
normalizable and reach the same correct RN-AdS asymptotic at spatial infinity. These extra solu- 
tions (which correspond to larger values of the grand canonical potential) lead to effective potentials 
that also vanish at the horizon and thus correspond to a non-zero conductivity at zero temperature. 
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I. INTRODUCTION 

The famous AdS/CFT correspondence flj allows to de- 
scribe confornial field theory in d-dimensional space-time 
by considering a. d-\- 1-dimensional super-gravity in anti- 
de Sitter space-time. This opens a number of opportuni- 
ties to look into non-perturbative quantum field theory 
at strong coupling. One of the recent interesting appli- 
cations of such a holography is constructing of a model 
of a superconductor. Usually in quantum field theory 
superconductors are well understood by the Bardeen- 
Cooper-Schrieffer, theory though there are indica- 
tions that for some systems the standard Fermi liquid 
theory cannot be a good approximation Therefore a 
holographic model for superconductors was suggested by 
HartnoU, Herzog and Horowitz Q . This model have been 
recently studied in a number of papers and some alter- 
native models of holographic superconductors were sug- 
gested fs^-fssj. These models contain a charged asymp- 
totically anti-de Sitter black hole which have non-trivial 
hairs at low temperatures. Until recent time, there were 
suggested various holographic models for the low tem- 
perature limit, while the dual description for the actual 
zero temperature ground state remained unknown. The 
very recent paper of Horowitz and Roberts solves 
this problem and find numerically the zero temperature 
holographic dual for superconductors. 

The system under consideration consists of the charged 
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scalar field coupled to a charged (3-|-l)-dimensional black 
hole, so that above some critical temperature, in the 
normal phase, the system is described by the Reissner- 
Nordstrom-anti-de Sitter black holes, while below the 
critical temperature, in the super-conducting phase, the 
black hole develops scalar hairs. Thus, the supercon- 
ductor is (2 -I- l)-dimensional, what might be realized 
for instance in graphene. In |35| . based on qualitative 
arguments, it has been shown that the effective poten- 
tial of the perturbation equation for the dynamic of the 
Maxwell field vanishes at the horizon, and, consequently, 
the conductivity never vanishes even at zero tempera- 
ture. Though some intuitive arguments were given in 
[35I about the behavior of conductivity in the suggested 
model, no calculations of conductivity were performed 
there, except for some estimations made for the low- 
frequency regime. Therefore our first aim here was to 
calculate conductivity for the Horowitz-Roberts model 
[ssj . When integrating the field equations, in addition 
to the ground state solution described in |35| . we have 
found a number of other solutions with the same leading 
AdS asymptotic at spatial infinity and obeying the same 
general form of anzats near the horizon. We have checked 
that the found here extra solutions, as that of JsS], have 
vanishing effective potential at the horizon, so that the 
conductivity will never be zero even at zero temperature. 
They correspond to configurations of the scalar field with 
higher energies at zero temperature. 

The paper is organized as follows. Sec II gives the ba- 
sics equations for the system of fields under consideration 
and scheme of construction of the numerical solution for 
a black hole with the scalar hair. Sec HI describes the 
spectrum of the obtained solutions which consists of the 
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ground state solution and solutions of the higher grand 
canonical potential. Sec IV is devoted to WKB calcula- 
tions of the conductivity for the zero-temperature super- 
conductor. 



II. CONSTRUCTION OF THE 
HOROWITZ-ROBERTS HOLOGRAPHIC DUAL 

The Lagrangian density for the system under consid- 
eration takes the form 

£ = i? + - -^F^'F^, - - iqA^j]' - Um (1) 

where V' is the scalar field, is the strength tensor 
of electromagnetic field, to, q are the scalar field's charge 
and mass and A is the vector-potential {F — dA). The 
cosmological constant is ~3/L^. The plane symmetric 
solution can written in a general form 

ds2 = -5(r)e-x('')dt2 + — + r2(dx2 + dy^) (2) 



A 



(r) dt, ip = ipii^) 



(3) 



We shall fix the gauge so that ip is real and measure all 
the quantities in units of the AdS radius, so that L ~ 1. 
The equations have the form: 

rJl-^A,^,f^,.-^^» ,4.) 

\g 2 rj 2g 



2 r 



2q^%l) 



-0 = 



(4b) 



1. The case = corresponds to a marginal oper- 
ator, A = 3, in the 2-1-1 superconductor with a 
nonzero expectation value. Following [sll we have 
used the ansatz 

<^ = r2+«, V = ^0 - ^^1^'^'+"^ (8) 
X = Xo-Xir^('+"\ .g-r2(l-.gir2(i+")). 

The coefficients in and g can be taken equal to 
unity. Substituting this into the field equations and 
equating the dominant terms for small r (with a > 
— 1), one has 
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+ 5a + 6 
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a'^ + 5a + 6 
4.{a + l) 



,Xo 



a + 2 



qe^ 



91 = 



+ 5a -I- 6 



(9) 



1/2 



2(2a2 + 7a + 5) V 2 

These formulas were obtained in [sBj. 

We solve the equations ^ numerically using the 
ansatz (jH)). We choose a in order to satisfy the 
condition using the shooting algorithm. 

2. In [s^ the ansatz for small r has been found for the 
case of TO^ < and g^ > —m?/& (A < 3): 

= 0or^(-ln(r))'/', V = 2 (- ln(r))'/' , (10) 
X = Xo + ln(- ln(r)), g = {2inn? /i)r^ ln(r), 

where 



, ,n rq-^(h'^^h^e^ , , 

X + + , = (4c) 



/3 



1-1¥>1. 



.'+(i-^)..+ ^-3.+ ^ = (4d) 

When we choose x = at infinity, the metric takes the 
standard AdS form at larger r. 

p y,(A) ^(3-A) 



where A = (3+V9 -f- 4to2)/2. In the boundary dual CFT, 
/i is the chemical potential, p is the charge density, and 
A is the scaling dimension of the operator dual to V'- We 
used 



^-(3-^) = 0. 



(6) 



The density of the grand canonical potential $7 of the 
state that corresponds to a given solution can be find by 
fitting the function g{r) at large r [3j| 



2 ^ 

g\T) = r H ho 



(7) 



Unfortunately, we were unable to construct a con- 
vergent procedure of integration in this case. We 
used the following method. We started the integra- 
tion of ([4]) from some point e which is very close to 
the horizon r = 0, substituting as an initial condi- 
tion the anzats (|10l) . Then we decreased the value 
of e and compare the results. We did not observe 
the convergence of the functions when decreasing 
e. Namely, as e approached zero the functions did 
not approach a certain limit, showing significantly 
different behavior. In the figure[T]we can see the de- 
pendence of the coefficient ij^'^^"^'^ on the parameter 
00 for various values of e. We checked that neither 
zeros of converge as e — and, therefore, we 

were not able to find the appropriate value of 0o for 
the solution that satisfies ^ . We have checked also 
that addition of sub-dominant terms to the anzats 
(|10p does not remedy the situation. 



Let us consider two cases: 



The TO = case is free from the above problem of 
absence of convergence and from here and on we shall 
consider only this case. 
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FIG. 1: Dependence of the coefficient ip^^ on cpo for = 
— 2 {X — 2) q — 2 for various staring points of integration 
e = 1/500 (cyan), e = 1/2000 (blue), e = 1/5000 (green), 
e = 1/20000 (red), e = 1/50000 (magenta). The smaller value 
of e is the closer zeros of ip^^~^^ are located. 



The equations (U) have the two-parametric symmetry 

M 

r ^ ar, t ^ -t, g -> a^g, (j) y4>, 6^e^. 
a b 

(11) 

One of these parameters allows us to fix x{t ^ °o) = 0. 
The other parameter re-scales the parameters of the so- 
lutions and can be chosen so that the chemical potential 
is unit. 

In order to satisfy these two conditions, after the solu- 
tion is found we choose 

5^e-x(oo)/2^ a = VA*- 

After this re-scaling the solution does not depend on 
the value of xo a-nd we can take Xo = 0- 

III. SPECTRUM OF THE SOLUTIONS 

In [3^ , a value of a satisfying ^ was found for m — 
as a function of q. It is interesting to note, that for 
each fixed q this value of a is not unique. At least for 
large values of q, there is a discrete spectrum of values 
of a. Each of these value of a (under the same fixed q) 
corresponds to a different solution of (|4]), that satisfies 
the condition The dependence q{a) for the first three 
solutions is shown on Fig. [2] We have checked that for all 
of the above three curves the solutions are normalizable 
and reach their AdS asymptotic at large distance. Near 
r = 0, all three solutions obey the same general anzats 
(|S]) though certainly with different values of a for each q. 
Although we have demonstrated only three solutions of 
the spectrum, it looks as if there is an infinite spectrum 
of solutions with increasing values of a for a fixed q. 

On the figure!?] we see how the coefficient -0(3-^) de- 
pends on a and on q: when a grows, the zeros of ?/)('^~-*') 



FIG. 2: Three lowest solutions for m — given by a as func- 
tions of q. 

become more and more dense in a, and when q grows 
the he zeros of ■0'^'^"^) become more spaced. Therefore 
different solutions (i.e. different lines q{ct)) lay closer to 
each other for smaller q and larger a making it difficult to 
distinguish numerically different nearby solutions. That 
is why the two upper curves do not continue on Fig. [5] to 
the region where they probably coincide or lay very close 
to each other: the numerical integration is not easy in 
that region as there are probably many other solutions 
nearby. We believe however that accurate numerical in- 
tegration could allow to complete at least a few upper 
curves until the minimal value of q. 

For larger a we observe that another singular point 
appears for r > (see Fig. [5]). 

The above found solutions correspond to lower energy 
states of the superconductor (see Fig. |3]). In order to see 
this, we find the density of the grand canonical potential 
for the corresponding states by fitting (O . Then, we can 
calculate the density of the free energy, using the formula 

F = f} + ^p. (12) 

We found that ft is larger for the state with larger a, 
but the free energy F appears to be lower (see caption 
for the figure 0. 

In order to check our numerical calculations of the 
thermodynamical potentials, we use the following rela- 
tion between the grand canonical potential and the free 
energy 

2n = -F. (13) 

This relation can be easily derived from (5.11) of [34] in 
the limit of zero temperature and magnetic field. We find 
that (|13p is satisfied up to the numerical precision. 

On the figure [2] one can see the three smallest values 
of a for which ^ is satisfied. The smallest a is the one 
found by Horowitz and Roberts. All the potentials are 
positive definite, vanish at the horizon {z = —00) and at 
the spatial infinity {z = 0). The potential for the lowest 
value of a has one peak. The effective potential for the 
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FIG. 3: The densities of the grand canonical potential (left figure) and the free energy (right figure) for the three lowest 
solutions for m = as functions of q. 




FIG. 4: Dependence of the coefficient ip^^ on a for q — & (left) and g = 30 (right). 



n-th higher value of a has n peaks (Fig. [6]). The larger 
value of a corresponds to the state with larger charge 
density plarger absolute value of the scalar hair ij}^'^'^ and 
lower density of the free energy. 



IV. CONDUCTIVITY BY THE WKB METHOD 

Assuming translational symmetry and stationary an- 
zats in time, the linearized perturbation of the vector 
potential satisfies the wave-like equation [s^ j 



A'L- 



A'M\^~—\e^- 1 A, = . 

(14) 

x/2 

Using a new radial variable dz = - — dr, at large r, dz — 
drjr^, and we choose the constant of integration so that 
z = — l/r. The horizon is located at z = — oo. Then (fT4)) 
has the wave-like form: 



(15) 



where the effective potential [35| 



(16) 



As was shown in |35| this effective potential always van- 
ishes at the horizon. Since we consider only solutions 
which satisfy ([6]), the potential also vanishes at the spa- 
tial infinity. 

In terms of non-rescaled functions the potential and 
the tortoise coordinate are given by 



V(r) 
dz 



a"- 
52 



dr 



a g 



-dr. 



(17) 



According to the Horowitz-Roberts interpretation, the 
holographic conductivity can be expressed in terms of the 
refiection coefficients j35| in the following way. In order 
to solve (|15p with the ingoing wave boundary conditions 
'At z — —00 we can extend the definition of the effective 
potential to positive z by setting y = for z > (the 
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g(r) 
0.5 - 




FIG. 5: Another singular point in the metric function g(r) 
at r > when a in the ansatz (|8]) is large (m = 0, g = 6, 
a = 2.8). 



boundary of the anti-de Sitter space (spatial infinity) is 
located at z = 0). 

Now an incoming wave from the right will be partly 
transmitted and partly reflected by the potential barrier. 
The transmitted wave is purely ingoing at the horizon 
and the reflected wave satisfies the scattering boundary 
conditions at z — oo. Thus the scattering boundary 
conditions for z > are 



= e~"^^ + Re" 
and at the event horizon 



-oo, 



(18) 



(19) 



where R and T are reflection and transmission coeffi- 
cients. Then one has 



A,{0) = 1 + R, A,^,iO) = -iujil-R). (20) 

3w: 
tivity is 



As shown in if = aI"'' + A^^ /r, and the conduc 



A 



(0) 



(21) 



Since A^^^ = -^^,^(0), so 



cr(w) = 



1 - R 
l + R 



(22) 



The above boundary conditions p8| . (|19p are noth- 
ing but the standard scattering boundary conditions for 
finding the S-matrix. The effective potential has the dis- 
tinctive form of the potential barrier, so that the WKB 
approach [sl] can be applied for finding R and a. Let 
us note, that as the wave energy (or frequency) w is real, 
the first order WKB values for R and T will be real [36j 
and 

+ ^ 1. (23) 
Next, we shall distinguish the two qualitatively different 
cases: first, when is much less then the maximum 
of the effective potential <C Vb, and second when 
is of the same order that the maximum of the potential 
cj^ ~ Vq and can be either greater or smaller than the 
maximum. Strictly speaking, we should have to consider 
also the third case when is much larger than the max- 
imum of the potential, but, as we shall see in most cases 
the reflection coefficient R decreases too quickly with w, 
so that a reaches its maximal value (unity) even at mod- 
erate UJ > Vq. 

For ~ Vq, we shall use the first order beyond the 
eikonal approximation WKB formula, developed by B. 
Schutz and C Will (see [36]) for scattering around black 
holes 



i?= (1 + e 



-2i7r(iy+(l/2)) 



UJ^^Vo, (24) 



where 



1 

^+2 



- Ao + A.I 



(25) 



Here Vq is the second derivative of the effective potential 
in its maximum, A2 and A3 are second and third WKB 
corrections which depend on up to 6th order derivatives 
of the effective potential at its maximum. 




FIG. 6: The effective potentials as functions of the corresponding tortoise coordinates for m = 0, g = 6 (A = 3) that for the 
four smallest values of a: 
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The above formula was extended up to the 6th WKB 
order in [s^] and applied to a number of problems of 
scattering around black holes (see for instance [111 and 
references therein). Mainly it was used for finding the 
so-called quasinormal modes of black holes, which im- 
ply special boundary conditions, so that ly becomes inte- 
ger in that case. For arbitrary v and each given uj the 
above WKB formula works for problems with the stan- 
dard scattering boundary conditions. We shall look for 
higher WKB orders in order to have the idea of possi- 
ble order of the error in the obtained results. Though 
the WKB series converges only asymptotically, in many 
cases, quite unexpectedly, WKB values have region of 



relative convergence in orders. 

The case of small frequencies is well described by the 
well-known formula 



T : 



(28) 



i.e. the at small frequency the transmission is exponen- 
tially suppressed. Here zi and Z2 are the turning points 



for which V{z) 
from (I23ll 



The reflection coefficient follows 



^2^^y^ (29) 

At small frequencies the reflection coefficient is close 
to 1, i.e. almost all energy is reflected by the poten- 
tial. Then R decreases with the increasing of w, and, for 
sufficiently large uj, usually seemingly larger than Vq or 
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FIG. 7: The conductivity found by the WKB formula for 
g = 10 (top, red), q = 12 (green), q — 14 (bottom, blue) as a 
function of frequency oj for m — 0. 



about it, the reflection coefRcient is close to zero. This 
means that according to (|^^ . the conductivity changes 
from zero at small frequencies until 1 at large frequencies. 
This kind of behavior we can see on Fig. [71 where the 
conductivity was obtained by using and the expression 
dUl) and the WKB formula (for w < 0.4) and IMD 
(for UJ > 0.4). There one can see that as approaches 
the peak of the potential barrier (which is located at 
r « 0.18), the accuracy of the formula ((29)) diminishes 
and ([M]) becomes a better approximation. Though a 
good confirmation of consistency of the both approxima- 
tions and ([M)) is the possibility of " smooth match- 
ing" of both data (see Fig. |8l |9]) if neglecting small in- 
termediate region of w, where both approximation have 
marginal accuracy (For Fig. 2, this intermediate region is 
0.3 < OJ < 0.4). In some range of parameters, such as the 
one shown on Fig. 3 (right), there is no such interme- 
diate region that should be neglected but, even better, 
both regimes ([^5]) and ([M)) overlap, giving almost the 
same values for some range of large values of q. 

Let us note here two important technical points. First 
is that when using formula ([2^ one needs the higher 
order derivatives of the effective potential which is un- 
known in analytical form, but is given only numerically. 
It would be a rough method to approximate the effec- 
tive potential by some interpolating analytical function 
and then to take derivatives of it: each derivative would 
bring additional numerical error to the calculations. In- 
stead we used the field equations (U) and have taken all 
necessary derivatives from ^ and by taking the corre- 
sponding derivatives of the wave and metric functions 4>, 
tfj, g, etc.. 

Another important moment is the accuracy of the used 
WKB technique. The existence of the "common region" 
where both formulas produce the same result says that 
for large values of q the WKB formulas work very well. 
The analysis of the higher order corrections indeed shows 
that for large q (and ~ Vb) the WKB series shows 



convergence in a few first orders: An example is q = 
10, w = 2.2, R = 0.778 for the first WKB order, R = 
0.722 for second WKB order, and R = 0.725 for the 
third order. This gives estimated error of less than one 
percent. There is no such good convergence for small 
values of q, therefore for ~ Vq we have used here the 
WKB formula of the first order for small q, and 3th order 
formula for large q. 

At small frequencies we have obtained a close, though 
not coinciding, numerical result to the formula (3.21) of 

M 



a= ( — ) , S = ^Wo + 1-1. 



(30) 



Thus for q = 1.6, we obt ained by WKB S w 4.8 (wq = 
0.45), what is close to y/Wo + 1 - 1 w 4.55, while for 
q = 1 WKB gives S « 3.0 (wo « 1.3) and V4V"o + 1 - 1 « 
3.97. The WKB correction to the as it can be seen 
from the above data, may be quite big and about 25 per 
cents. 

The formula pop used for its derivation rough match- 
ing of the left and right dominant asymptotics. Thus it 
is expected to be less accurate than the WKB method 
we are using here. 



V. CONCLUSIONS 

We have found the WKB values of conductivity for the 
Horowitz-Roberts model of the zero-temperature super- 
conductor for = case. Dependence of conductivity 
on parameters of the theory such as the charge density p 
and the frequency w is investigated. WKB data for con- 
ductivity confirms the qualitative arguments that a does 
not reach zero even at zero temperature, in agreement 
with [3^. By the WKB calculations we have confirmed 
the analytic relation derived in [35] for the w-dependence 
of a at small frequencies and calculated the pre-factor 
for this relation for various q. The used here third order 
WKB formula which has very good accuracy for large val- 
ues of q (and moderate w), showing convergence in orders 
with an estimated error of around fractions of one per- 
cent. In addition, we have found the set of other solutions 
which describe the superconductor at zero temperature 
in the states with higher grand canonical potential and 
lower energy. 

Our paper may be improved in a number of ways. First 
of all, the conductivity values could be obtained with bet- 
ter accuracy, if one uses the numerical shooting, which is 
known to work well for asymptotically AdS space-times 
[43 |. The conductivity of higher grand potential states 
with m — cannot be obtained by WKB formula we 
used, because the effective potentials for higher states 
have a number of local maximums. Thus accurate shoot- 
ing approach would allow also for complete analysis of 
conductivities of these states. Finally, the case of non- 
vanishing mass of the scalar field m probably requires 
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FIG. 8: The conductivity for m = 0, g = 1 as a function of frequency lj (5 ~ 3.0). On the left figure the sohd line corresponds 
to the interpolation between the two WKB approximations (I24p and (|29[1 . On the right figure the solid line corresponds to the 
fit of the numerical data for small values of cj. 
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FIG. 9: The conductivity for m = 0, q = 1.6 as a function of frequency uj (5 ~ 4.8). On the left figure the solid line corresponds 
to the interpolation between the two WKB approximations (|24p and (|29p. On the right figure the solid line corresponds to the 
fit of the numerical data for small values of uj. 



FIG. 10: WKB conductivity for m — and for various values q uj — 1/2 (left) and uj — 5 (right). As q grows, the maximum of 
the effective potential grows, so that a fixed uj moves down from the peak, approaching the regime iv'^ <^Vo (red dots). 
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some other and more sophisticated procedure of integra- 
tion or a different anzats near the horizon from that sug- 
gested in [ssj ]. 
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